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Nonlinear evolution of a trans- Alfvenic shock wave (TASW), at which the flow velocity passes 
over the Alfven velocity, is computed in a magnetohydrodynamic approximation. The analytical 
theory suggests that an inflnitesimal perturbation of a TASW results in its disintegration, i.e., finite 
variation of the flow, or transformation into some other unsteady conflguration. In the present 
paper, this result is confirmed by numerical simulations. It is shown that the disintegration time is 
close to its minimum value equal to the shock thickness divided by a relative velocity of the emerging 
secondary structures. The secondary TASW that appears after the disintegration is again unstable 
with respect to disintegration. When the perturbation has a cyclic nature, the TASW undergoes 
oscillatory disintegration, during which it repeatedly transforms into another TASW. This process 
manifests itself as a train of shock and rarefaction waves, which consecutively emerge at one edge 
of the train and merge at the other edge. 



I. INTRODUCTION 

It has long been believed that trans-Alfvenic shock 
waves (TASWs), at vifhich the flow velocity passes over 
the Alfven velocity, cannot exist in the real world. Since 
a stationary trans-Alfvenic shock transition was obtained 
in a numerical simulation [0 , this conventional view point 
was replaced by an opposite view point. The overall 
claim was that there is no principal difference between 
TASWs and fast and slow shocks, at which the flow is 
super- and sub- Alfvenic, respectively. At the same time, 
the contradiction inherent in a stationary TASW, which 
follows from an analytical theory, was not lifted. To rec- 
oncile this contradiction, it was suggested that a TASW 
exists in an unsteady state in which it is repeatedly de- 
stroyed and recovered In the present paper, we show 
by way of magnetohydrodynamic (MHD) simulation that 
the evolution of a TASW may have the form of oscil- 
latory disintegration, i.e., reversible transformation into 
another TASW. 

The disintegration of an arbitrary hydrodynamic dis- 
continuity was considered for the flrst time by Kotchine 
1^. After that, Bethe [Q studied the disintegration of 
shock waves. In the absence of a magnetic field, the shock 
may disintegrate only in a medium with anomalous ther- 
modynamic properties. The magnetic field enlarges the 
number of possible discontinuous structures thus giving 
additional degrees of freedom for the disintegration. The 
disintegration configurations of arbitrary MHD disconti- 
nuities were obtained in Refs. Q. Furthermore, it has 
been shown that trans-Alfvenic shock transitions can be 
realized also through a set of several discontinuities , 
in contrast with fast and slow transitions. However, this 
fact on its own does not assure that the shock disinte- 
grates. 



The important feature that predetermines the disinte- 
gration of TASWs is their nonevolutionarity. The prob- 
lem of evolutionarity was initially formulated for the 
fronts of combustion |^] and hydrodynamic discontinu- 
ities . Evolutionarity is a property of a discontinuous 
flow to evolve in such a way that the flow variation re- 
mains small under the action of a small perturbation. 
This is not the case for a nonevolutionary discontinuity. 
At such a discontinuity, the system of boundary condi- 
tions, which follow from the conservation laws, does not 
have the unique solution for the amplitudes of outgoing 
waves generated by given incident waves. From a math- 
ematical view point this means that the number of un- 
known parameters (the amplitudes of the outgoing waves 
and the discontinuity displacement) is incompatible with 
the number of independent equations. Since a physical 
problem must have the unique solution, the assumption 
that the perturbation of a nonevolutionary discontinu- 
ity is inflnitesimal leads to a contradiction. In fact, the 
infinitesimal perturbation results in disintegration, i.e., 
finite variation of the initial flow, or transformation into 
some other unsteady configuration. 

The evolutionarity requirement gives additional re- 
strictions on the flow parameters at a shock, compared 
to the condition of the entropy increase. The restrictions 
appear because the direction of wave propagation (to- 
ward the discontinuity surface or away from it), and thus 
the number of the outgoing waves, depends on the flow 
velocity. If the velocity is large enough then the given 
wave may be carried down by the flow. Therefore, at an 
evolutionary discontinuity, the flow velocity must be such 
that it provides the compatibility of the boundary equa- 
tions. This form of evolutionarity condition was applied 
to MHD shock waves in Refs. HJl^. As a result, the fast 
and slow shocks are evolutionary, while the TASWs are 
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nonevolutionary. 

This classical picture was challenged when Wu |jl|] ob- 
tained a stationary TASW in a numerical simulation. 
The existence of a stationary numerical solution does not 
mean of course that the shock is stable with respect to 
disintegration or transition into another unsteady flow. 
Wu ||T|] demonstrated that a TASW, which is subfast 
upstream and subslow downstream, disintegrates under 
the action of a small Alfven perturbation with a large 
enough characteristic time. Nevertheless, this numerical 
result was interpreted as being in a contradiction with 
the principle of evolutionarity and stimulated the efforts 
to modify or even disprove this principle. 

It was suggested that the free parameters that describe 
a nonunique structure of a TASW ||l^ or the amplitudes 
of strongly damping dissipative waves [ p^ should be in- 
cluded in the number of unknown parameters when solv- 
ing the problem of evolutionarity. This would make the 
TASW evolutionary. In both cases, however, the per- 
turbation is confined within the shock transition layer. 
Consequently, it does not enter into the boundary con- 
ditions, which relate the quantities far enough from the 
transition layer, and thus it does not contribute to the 



evolutionarity |14|. 

Wu |lj] also argued that the TASW whose nonevolu- 
tionarity is based on separation of Alfven perturbations 
from the remaining perturbations | [Tc| ] becomes evolution- 
ary in the case of a nonplanar shock structure because 
in this case the separation formally does not take place. 
However, as shown by Markovskii iQ, the coupling of 
the small-amplitude Alfven modes with a low enough fre- 
quency to the remaining modes is weak (unless the shock 
is of the type close to one of the degenerate types, Alfven 
discontinuity or switch shocks). Therefore the coupling 
becomes essential only when the small perturbation gen- 
erates large variation of the flow, which is the same result 
as predicted by the principle of evolutionarity. 

There is one more finding that favors the nonexis- 
tence of stationary TASWs. As discussed by Kantrowitz 
and Petschek ||l^, the TASWs are isolated solutions of 
Rankine-Hugoniot problem, which do not have neighbor- 
ing solutions corresponding to small deviations of bound- 
ary conditions. Wu and Kennel introduced a new 
class of trans-Alfvenic shock-like structures with non- 
coplanar boundary states. The thickness of such a struc- 
ture increases in the course of time, and it eventually 
evolves to a large-amplitude Alfven wave. It was thus 
shown that neighboring to a TASW are time-dependent 
configurations, which are not solutions of the Rankine- 
Hugoniot problem. In addition, Falle and Komissarov 
I prt recently considered stationary TASWs of all possi- 
ble types and showed that the shocks disintegrate if the 
boundary values deviate from their initial values. 

Strictly speaking, a TASW, in contrast with fast and 
slow shocks, becomes a time-dependent shock-like struc- 
ture once it is perturbed by a small-amplitude Alfven 
wave because the Alfven wave violates the coplanarity 
condition. This fact, already on its own, means that the 



TASW becomes unsteady under the action of the small 
perturbation. However, the scenario for its evolution de- 
pends on the initial configuration and on the nature of 
the perturbation. After the disintegration, the magnetic 
field reversal given at the initial nonevolutionary shock 
may be taken either by a secondary TASW or by an 
Alfven discontinuity. Both structures are nonevolution- 
ary H;|l^. Therefore single disintegration does not lift 
the contradiction inherent in a TASW. The main ques- 
tion that we solve in this paper is what happens to the 
post-disintegration nonevolutionary configuration. We 
show that the secondary TASW is again unstable with re- 
spect to disintegration and that the evolution of a TASW 
may have the form of oscillatory disintegration. In Sec. 
II, we describe the simulation method. In Sec. HI, we 
discuss the results of the calculations. Our conclusions 
are presented in Sec. IV. 



II. NUMERICAL METHOD 

We take the MHD equations in the following form 
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Here the subscript "_L" denotes the vector component 
perpendicular to the x axis, Bx = const, magnetic diffu- 
sivity I'm and viscosity 77 are put constant and equal to 
0.1 in all calculations, and we use the units such that the 
factor 47r does not appear. The initial distribution of the 
MHD quantities is given by the following formulas 



hiPl + Pi) ^ UpI ^ Pi)tanh(a;/L), 



(2a) 



Vx = ^{vx] + Vxi) - ^{vx] - Vxi)ta.nh{x/L), (2b) 



P= kiP^ +Pl) - Up] --Pi)tanh(a:/L), (2c) 
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By = BrCos{e), B^=BrSm{9), (2d) 

Vy = VrC0s{6), Vz = VrSin{6), (2e) 

Br = i(| S±T I + I I) 

-1{\B^^\-\B^^ |)tanh(x/L), (2f) 

= 5(1 W_LT I + I I) 

-i(| t^xT I - I t'ii |)tanh(a;/L), (2g) 

6* = f (1 + tanh(a;/L)), (2h) 



where the subscripts "t" and "|" denote the quantities 
in the asymptotic upstream and downstream regions, re- 
spectively. 

After the configuration relaxes to a steady state, it is 
perturbed by an Alfven wave specified by the expression 



B,^^dBzi^l + tanhi^^-^jj, (3a) 

Vz = -Bzl^p. (3b) 

This wave moves downstream. The configuration is set 
by putting Bj. = 0.89, L = 1.45, and 

By-^ = 0.93, Byi = -0.8, (4a) 

S,T-t;,t=0., Bzi^Vzi=Q., (4b) 

u^l = 1., = 0.55042, (4c) 

Vy^ = 1.04494, Vyi = -0.49476, (4d) 

P\ = 1., Pi - 1.81681, (4e) 

Pt = 0.00116, Vi = 0.56319. (4f) 



This corresponds to a II ^ III shock, for which V+| > 
Vx\ > VaxI and VaxI > Vxi > V-i, where V+ and are 
the fast and slow magnetosonic velocities. 

We solve Eq. (|l]) using a uniform grid and an ex- 
plicit conservative Lax-Wendroff finite-difference scheme 
with physical dissipation [|l8|. The time step is limited 
by the Courant-Friedrichs-Lewy (CFL) condition and by 
the dissipation timescale. The boundary values are ob- 
tained by hyperbolic interpolation. The numerical inter- 
val —50 < a; < +300 is covered by 2600 grid points. The 
interval is chosen in such a way that no large-amplitude 
wave reaches the boundaries during the computation 
time. Small-amplitude waves pass through the bound- 
aries without any detectable reflection which could affect 
the flow inside the simulation region. We have tested our 
code for a smaller mesh and a corresponding time step 
determined by the CFL condition as well as for the same 
mesh and a time step smaller than that determined by 
the CFL condition. The test showed that there is no con- 
siderable dependence of our results on the mesh size and 
time step. 



III. RESULTS OF SIMULATIONS 

Equation @) does not exactly describe the shock struc- 
ture. Therefore the flow undergoes time variations until 
it adjusts to a stationary shock transition. The result- 
ing boundary values differ slightly from those given by 
Eq. (^) but the difference is less than 1%. The con- 
servation laws for these new values are fulflUed with the 
precision less than 0.1%. The stationary configuration is 
then perturbed by a small-amplitude Alfven wave with 
I = -L, xo = -40, and SBz = 0.025 or SBz = -0.025 
(Fig. |l]). Note that SBz is about 50 times smaller than 
I B| I . Although in the case of an upstream incident 
wave the perturbation of Bz and Vz (not shown) is carried 
to the downstream region, the boundary conditions for 
the Alfven waves are incompatible. Therefore the given 
Alfven perturbation pumps Bz and Vz into the shock or 
out of the shock, depending on the sign of Bz inside the 
transition layer. Since Bz inside the transition layer is 
nonzero, the shock behaves in different ways under the 
action of the perturbations with positive and negative 
SBz- If the shock and the perturbation carry Bz of the 
same sign, the shock disintegrates into a II — > III shock 
of a smaller amplitude, a large-amplitude slow shock, 
and some other structures of a much smaller amplitude 
(Fig. |a,b). 

If the shock and the perturbation carry Bz of oppo- 
site signs, the situation is somewhat peculiar. The main 
secondary structures are a TASW and a slow rarefaction 
(Fig. ||a,b). However, these structures do not become 
separated. The reason is that the secondary TASW is of 
a so- called II ^ IV = III type |l^. This means that the 
downstream velocity at the shock is exactly equal to the 
slow magnetosonic velocity. Therefore there is no disinte- 
gration in the usual sense but the configuration becomes 
unsteady because the right boundary of the slow rarefac- 
tion moves away from the TASW, while the left boundary 
remains attached to the shock. Note that the rarefaction 
wave is attached to the TASW not at the density peak 
but somewhere to the right of the peak. This is related 
to the fact that the density profile of a II ^ IV shock 
has a maximum (see, e.g., Rcf. pT| ), in contrast with the 
monotonic profile of a II ^ III shock. 

From the moment when the Alfven wave with SBz > 
arrives to the shock, the disintegration starts almost im- 
mediately, in contrast with the result of Wu [0 . The rea- 
son is that the disintegration time depends on the shock 
type and on its initial state. This can be understood as 
follows. The important characteristic of a TASW, intro- 
duced by Kennel et al. |jl^ , is the integral of Bz over the 
transition layer, 

Iz^ J Bzdx. (5) 

This integral fixes the nonunique structure of a TASW. 
For a II ^ III shock, the quantity Iz takes two distinct 
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values, Izo and —Izo, and, for a I ^ III or II IV shock, 
it falls into the interval —Izo < Iz < IzO- The quantity Izo 
depends on the boundary values, and it tends to infinity 
when the shock approaches an Alfven discontinuity or a 
switch shock, which is intermediate between evolutionary 
and nonevolutionary shocks. This result was obtained 
for almost parallel small- amplitude shocks, but one may 
expect that it remains qualitatively valid in the general 
case. 

When an Alfven wave is incident on a TASW, it 
changes Iz ■ If we start from a planar I ^ III or II IV 
shock {Iz — 0), as in the case studied by Wu Q, the 
quantity | Iz \ first has to reach the value Izo- Only after 
that it falls into the forbidden region, and the disinte- 
gration starts. In the case of a II ^ III shock, there is 
a different situation. Since Iz takes only distinct val- 
ues Izo and —IzQ, the disintegration starts immediately, 
and the disintegration time is close to its minimum value 
L/V, approximately equal to 30 in our case, where V is 
a relative velocity of the secondary discontinuities. 

Let us now follow the further evolution of the post- 
disintegration configuration under the action of a small 
perturbation. Our main conclusion is that the secondary 
TASW is again unstable with respect to disintegration. 
At the same time, the way of evolution depends on a 
form of the perturbation. We first discuss the case where 
the perturbation of the secondary TASW is such that 
Iz continues to increase or decrease, in particular where 
the perturbation is equal to its initial positive (Fig. |^b,c) 
or negative (Fig. |pD,c) value. If the perturbation of Bz 
is positive then the shock spreads in space, with all the 
jumps, except for ABy and Avy, decreasing in time. It 
thus approaches a large-amplitude Alfven wave. If the 
perturbation is negative, the shock first passes through 
the state in which Iz — 0. This is not in a contradiction 
with the analytical theory, because a II ^ IV shock may 
have a planar structure, in contrast with a II — s- III shock 
pof . When | Iz \ reaches a critical value, the shock disin- 
tegrates (Fig. ^), and after that it spreads in space ap- 
proaching a large-amplitude Alfven wave (Fig. ^). The 
precursor of the disintegration is the peak in By curve in 
Fig. §3. 

We now turn to a cyclic perturbation. We impose the 
perturbation described by Eq. (|^) in such a way that 
Bz changes sign at a; = —40 and the Alfven wave now 
carries the perturbation of the same amplitude but oppo- 
site sign. After the first disintegration starts (at t = 20 
for SBz > and at t = 670 for 5Bz < 0), the opposite 
sign perturbation arrives to the TASW each 150 units 
of time. The resulting configuration is such that the in- 
crease of I /z I is repeatedly replaced by its decrease, and 
the shock undergoes oscillatory disintegration. The dis- 
integration configurations after several cycles are shown 
in Figs. H and ^. As can be seen from the figures, the 
configuration emits a sequence of contact discontinuities. 
The contact discontinuities move with the fiow velocity, 
which is approximately equal to that given by Eq. (^). 
The corresponding time interval between the discontinu- 



ities is equal to 150. 

Downstream of the TASW, there is a wave train, which 
consists of slow shock and rarefaction waves. These struc- 
tures are not standing in the fiow. They consecutively 
emerge at the left edge of the train and merge at the 
right edge. The merging is seen in Fig. ^ at a; w 150. 
We note that, in the case of a negative initial perturba- 
tion, the transition through the state with = is not 
necessary for the oscillatory disintegration to occur. If 
the perturbation changes sign for the first time before 
Iz becomes negative, the disintegration configuration is 
similar to that shown in Fig. ^, except for the sign of Iz 
inside the shock. 

Finally, the shock comes to a steady state only in a 
degenerate case where the perturbation of the secondary 
TASW exactly compensates the nonzero value of Bz and 
Vz outside of the transition layer. We emphasize that 
in all but the degenerate cases the small Alfven pertur- 
bation makes the TASW unsteady, in contrast with fast 
and slow shocks. However, there remains a question. For- 
mally, the initial TASW becomes a time-dependent struc- 
ture, much like the secondary TASW, since the Alfven 
perturbation arrives to the initial shock. The question 
is why the initial TASW disintegrates when Iz increases 
monotonically, while the secondary TASW does not. To 
answer this question, we first mention that the secondary 
TASW is more close to a finite-amplitude Alfven wave 
than the initial shock. Alfven waves, as well as switch 
shocks, are singular structures. As shown by Kennel et 
al. , the quantity dizo /dq tends to infinity as the shock 
approaches these singular structures. Here q = By^/Byi 
characterizes the jumps of the boundary values at the 
shock with a given Iz, and Izo{q) is an allowed curve in 
which a II — > III shock has a stationary structure. 

Assume now that the initial shock is in the state 
Iz — Izoilo)- A small Alfven perturbation changes Iz- 
For the shock to remain in the curve Izo{q), a change 
of q is required. In the general case, the variation of 
Iz is comparable with the variation of q, and thus with 
the jumps of the boundary values at the TASW. In this 
case, the evolution has the form of disintegration. By 
contrast, if the shock is close to the singular structure, 
the given variation of Iz requires a small variation of q, 
and the jumps of the boundary values are adjusted to 
Izoil) in a diffusion- like manner. It should be mentioned 
that the curves Izo{q) were obtained by Kennel et al. 
JlOt for small-amplitude shocks propagating almost par- 
allel to the magnetic field. Nevertheless, we speculate 
that, in our simulation, the initial TASW has a small 
enough dizo/dq to disintegrate, while for the secondary 
TASW the quantity dlzo/dq is large enough to dim the 
disintegration. Such an explanation does not imply that 
a TASW cannot disintegrate more than one time in prin- 
ciple. Furthermore, in our simulation run with a positive 
constant SBz, there is an evidence for a possible second 
disintegration at t — 380. However, the second disinte- 
gration is too faint to contend that it indeed takes place. 
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IV. CONCLUSIONS 
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We have performed a numerical simulation of a trans- 
Alfvenic shock wave. The shock that we have considered 
is of a II ^ III type, i.e., it is subfast upstream and super- 
slow downstream. We have shown that the shock disin- 
tegrates under the action of a small Alfven perturbation. 
The resulting configuration includes a secondary TASW, 
a large-amplitude slow shock or rarefaction wave, and 
other small-amplitude structures. We have also demon- 
strated that the secondary TASW is again unstable with 
respect to disintegration. When the perturbation has a 
cyclic nature, the shock undergoes an oscillatory disin- 
tegration. This result is in a qualitative agreement with 
our previous finding . This process shows up as a train 
of slow shock and rarefaction waves, which consecutively 
emerge at one edge of the train and merge at the other 
edge. At the same time, the disintegration configuration 
of a small-arnplitude almost parallel TASW discussed by 
Markovskii |^| includes alternating TASWs and Alfven 
discontinuities rather than alternating TASWs. This dis- 
crepancy is explained by the fact that, in the approxima- 
tion used in Ref . |^ , the diS'erence between the secondary 
TASW and the Alfven discontinuity manifests itself in 
higher orders. 

In contrast with the results of Wu , the disintegra- 
tion starts almost immediately after the Alfven perturba- 
tion arrives to the initial shock. The characteristic time 
of this process is equal to that required for the secondary 
structures to become separated. The reason for this can 
be seen as follows. TASWs have a nonunique structure. 
A II ^ IV shock transition studied by Wu as well 
as a I — > III transition, allows a continuous family of in- 
tegral curves, while the II — > III shock has two distinct 
integral curves. For given boundary values, each integral 
curve is fixed by the definite parameter. The incident 
Alfven wave changes the parameter and thus the shock 
structure. In the case of a I — > III or II IV shock, some 
time passes until the parameter falls into a forbidden re- 
gion, and only after that the shock disintegrates. In the 
case of a II — > III shock, its structure immediately be- 
comes inconsistent with the boundary values under the 
action of the Alfven wave, which initiates the disintegra- 
tion. 

Thus, our simulations confirm that a TASW becomes 
unsteady when it is perturbed by a small-amplitude inci- 
dent wave. Furthermore, an almost vanishing perturba- 
tion results in considerable dynamics at relatively small 
timescales. The scenario for the shock evolution depends 
on its initial state and on the nature of the perturba- 
tion. In particular, the evolution may have the form of 
oscillatory disintegration in which the shock repeatedly 
transforms into another TASW. 
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FIG. 1. Inintial distribution of the density and magnetic 

field at the TASW perturbed by a small-amplitude Alfvon 
wave with positive (solid line) and negative (dashed line) value 
of SB^. 
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FIG. 2. Disintegration configuration for a constant positive 
perturbation at t = 100 {a),t = 260 (b), and t = 700 (c). 
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FIG. 3. Disintegration configuration for a constant nega- 
tive perturbation at t = 100 {&), t = 700 (b), and t = 1100 
(c). 
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FIG. 4. Configuration after several cycles of duration 150 
time units at t = 975 {&), t = 1125 (b), and t = 1275 (c). 
Oscillatory disintegration is started by a positive initial per- 
turbation. The perturbation changes sign for the first time at 
t = 150. 
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FIG. 5. Configuration after several cycles of duration 150 
time units at t = 1300 (a), t = 1450 (b), and t = 1600 
(c). Oscillatory disintegration is started by a negative initial 
perturbation. The perturbation changes sign for the first time 
at t = 800. 
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